Large quantum nonlinear dynamic susceptibility of single-molecule magnets 
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The nonlinear dynamical response of Mni2 single-molecule magnets is experimentally found to 
be very large, quite insensitive to the spin- lattice coupling constant, and displaying peaks reversed 
with respect to classical superparamagnets. It is shown that these features are caused by the strong 
field dependence of the relaxation rate due to the detuning of energy levels between which tun- 
neling takes place. The nonlinear susceptibility technique, previously overlooked, is thus proposed 
as a privileged probe to ascertain the occurrence of quantum effects in mesoscopic magnetic systems. 
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The possibility to observe quantum phenomena, such 
as tunneling or even quantum coherence, occurring in 
mesoscopic magnetic or superconducting materials is a 
subject of considerable interest Q. These systems are 
also important for the study of the fuzzy border be- 
tween the classical and quantum worlds @. As the size 
increases, tunneling probabilities decrease exponentially 
and the system becomes more entangled with its envi- 
ronment (phonons, electromagnetic fields, nuclear spins, 
electrons, etc). This interaction plays a crucial role as it 
enables the system to approach equilibrium but also 
causes decoherence of its quantum dynamics |4j . Besides 
its fundamental interest, decoherence sets the limits for 
application of the quantum behavior of mesoscopic sys- 
tems to, e.g., quantum computing. 

Molecular crystals are ideal for these studies . They 
are ordered lattices of identical molecules, each contain- 
ing a magnetic cluster that is isolated from its neigh- 
bors by a shell of organic ligands. Because of strong 
intramolecular super-exchange interactions, each cluster 
behaves at low temperatures as a superparamagnet with 
a large spin S that flips as a whole. The best studied of 
these is the Mni2 acetate, for which S — 10. Its strong 
magnetic anisotropy produces magnetic relaxation and 
hysteresis at low temperatures which make it a po- 
tential candidate for information storage at the molecular 
level. The spin reversal takes place via quantum tunnel- 
ing through the anisotropy energy barrier whenever the 
magnetic field brings into resonance magnetic states lo- 
cated at both sides of the anisotropy barrier [t| • 

For classical superparamagnets, Garcfa-Palacios and 
Svedlindh §| predicted that the nonlinear dynamic sus- 
ceptibility can be very large and, in contrast to the linear 
susceptibility, quite sensitive to the spin-bath coupling- 
strength (damping parameter). In this work, we apply 
this technique, unexploited in this field, to investigate en- 
vironmental effects in molecular clusters. We have stud- 
ied two compounds, both having the same Mni2 magnetic 
core, but different ligand molecules. This enables us to 



investigate how the spin dynamics is affected by different 
couplings to the phonon bath. 

Our results show that the nonlinear response of these 
single-molecule magnets is dominated by a very large 
contribution, hitherto unforeseen, whose sign is opposite 
to the classical one. This contribution is interpreted in 
terms of the extreme sensitivity of quantum tunneling 
probabilities to the detuning of the magnetic levels by the 
external field. Therefore the dynamical nonlinear suscep- 
tibility permits to ascertain unambiguously if tunneling 
plays a role in the relaxation of magnetic molecular clus- 
ters and nanoparticles. In addition to its fundamental 
interest, this phenomenon might be relevant in situations 
where nonlinear effects become important (e.g., magneto- 
optical applications of these systems). 

The preparation of the Mni2 acetate, 
[Mni2(CH 3 COO)i 6 (H20)40i2]-2CH3COOH-4H20, 
was based on the original method of Lis Q . The Mni2 2- 
Cl benzoate, [Mn^^ClPhCOO^e^O^Oial^Chj ■ 
5H2O, hereafter abbreviated to M1112CI, was prepared 
following the method described in Ref. ^(J- XRD 
patterns of powders of both materials show peaks 
corresponding to the crystal structures described in 
the literature 0, 0|. In order to avoid uncertainties 
concerning the orientation of the crystals, we measured 
powdered samples embedded in a two-component glue. 

The nonlinear susceptibility is usually obtained (see, 
e.g., Ref. from the third-harmonic of the response 
to an ac magnetic field ~ Ahe luJt . Alternatively |l2j, one 
can consider the first-harmonic response x( w ) but in the 
presence of weak superimposed dc fields H, from which 
one can get a quantity that posseses similar dependences 
on the relevant parameters (frequency, damping, etc.). 
Specifically, the (first) nonlinear susceptibility \3 was ob- 
tained by fitting the real and imaginary components of 
the field-dependent ac susceptibility to a polynomial 



X(w, H) = xi H + 3x 3 M# 2 + 5 X5 MH 4 + 



(1) 



where xi is the custom linear susceptibility. The "expan- 
sion coefficients" Xi;X3;--- are H- independent and the 
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FIG. 1: Dynamical susceptibility of M1112CI, normalized by 
its zero-field value, measured at T = 4.5 K as a function of 
the dc magnetic field. Full symbols, real part; open sym- 
bols, imaginary part. The frequencies shown correspond to 
lot = 0.015, 0.4, and 265. The solid lines represent polyno- 
mial fits (of order 8) from which the non-linear susceptibility 
is obtained [see Eq. Jy]. The dotted line is the parabolic term 
Xi + 3x3^7 2 that dominates the low-field behavior. 



factors 1, 3, 5, . . . arise when raising to n in Xn(H + Ah) n 
and retaining terms linear in Ah. Below we will show 
that the classical X3 obtained in the above way is equiv- 
alent to the third-harmonic susceptibility. 

All magnetic measurements were performed using a 
commercial SQUID magnetometer. Measurements were 
performed in a broad range of temperatures 1.8 K < T < 
60 K and frequencies 0.025 Hz < w/2ir < 1.5 kHz. We 
took especial care to ensure that the amplitude of the ac 
field (Ah = 4.5 0c) was sufficiently small not to induce 
any additional nonlinearity (i.e., Xi>X3> ■ ••■> _ are indepen- 
dent of Ah). Previous results on x(<J, H) [14J contain the 
dependences of the linear and nonlinear susceptibilities 
mixed, whereas the above method enables one to sepa- 
rate and analyze the genuine nonlinear contributions. An 
example of its application is shown in Fig. ^ 

Linear ac susceptibility data measured as a function of 
to at zero field are shown in Fig. [21 The results are well 
described by a Debye law 



X = Xs + 



XT ~ XS 

1 + i LOT ' 



(2) 



where xt and xs are the isothermal (equilibrium) and 
adiabatic limits of xi- The samples have therefore well- 
defined relaxation times r, which follow an Arrhcnius 
law t — To exp (U /k&T). The activation energy U is the 
same for both compounds (same slopes in the inset of 
Fig. EJ, i.e., tunneling takes place via the same excited 
levels, corresponding approximately to the tunnel-split 
m = ±4 states. In contrast, the pre-exponential factor tq 
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FIG. 2: Real (circles) and imaginary (triangles) parts of the 
reduced linear susceptibility (xi — Xs)/(xt — Xs), measured 
between T = 3.7 K and 6.6 K in the frequency range 2.5 x 
10" 2 Hz < uj/2ti < 1.5 x 10 3 Hz. The lines are calculated 
using Debye law [Eq. The inset shows the temperature 

dependence of the relaxation time r. Full and open symbols 
correspond to Mni2 acetate and Mni2Cl, respectively. 



of Mni2Cl is a factor two shorter than tq of Mni2 acetate. 
To is related to the lifetime of excited states, which decay 
by emitting phonons to the lattice. Therefore, the spin- 
lattice interaction, or damping, appears to be stronger 
for Mni2Cl, in agreement with previous results This 
shows that these two materials are suited to study how 
a magnetic cluster behaves in different environments. 

Before presenting the results for the nonlinear suscep- 
tibility for these systems, let us consider the theoretical 
behavior of this quantity in the classical limit. We have 
calculated X3 from the exact (continued-fraction) solu- 
tion of the corresponding Fokker-Planck equation 
To mimic the experiments: (i) the calculation is done 
for spins with randomly distributed anisotropy axes and 
(ii) the first harmonic of the response is computed in 
weak dc fields, and the quadratic coefficient extracted 
[see Eq. JIJ]. The results are displayed in Fig.|31for vari- 
ous values of the Landau-Lifshitz damping parameter A. 
It is seen that the ui dependence of the so-obtained X3, as 
well as its damping sensitivity, are qualitatively similar to 
those of the third harmonic X3 studied in Ref . || . In the 
classical case A determines the time available for the spin 
to precess around the anisotropy axis, before decaying to 
the bottom of the potential wells. 

The experimental X3 data of both compounds are also 
shown in Fig. (lower panel). They agree, when lot <C 1 
(sufficiently high temperatures or low frequencies), with 
the expected equilibrium behavior X3T l/T 3 . When 
~ 1) X3 becomes strongly dependent on to. It ex- 
hibits a large maximum X3 — 6OX3T near lot = 1, which 
depends weakly on the damping (i.e., on To). However, 
even more striking is the fact that the experimental and 
classically calculated peaks of X3( w ) have opposite signs! 

We next try to explain the origin of this large and 
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FIG. 3: Upper panel: nonlinear susceptibility of classical 
spins, calculated by the method of Ref. [g, at ksT '/U = 5/70 
for three damping coefficients A. Full lines, real part; dotted 
lines, imaginary part. Lower panel: nonlinear susceptibility 
of Mni2 acetate (full symbols) and Mni2Cl (open symbols) 
measured at T = 5 K. Dots and triangles represent the real 
and imaginary parts. We have multiplied to by the zero field 
t = 1.31 x 10~ 2 s and 5.55 x 10~ 3 s, for Mni2 acetate and 
Mni2Cl, respectively. The lines are quantum calculations per- 
formed as in Ref. \L7\: full line, real part; dotted line, imag- 
inary part. The inset shows the relaxation time of a single 
crystal of Mni2 acetate measured with the field applied along 
the anisotropy axis (•) and the classical result (full line). 



(classically) unexpected nonlinear response. We observe 
that Eq. J2J a l so describes the susceptibility of Mni 2 in 
moderate fields. Therefore, X3 can be estimated by differ- 
entiating Debye's expression with respect to H, getting 



X3 



X3T 



1UJXT 



6(l + iur) 2 dH 2 



(3) 



where the subscript indicates evaluation at H = 0. The 
second, dominant, term is directly related to the field 
derivative of the relaxation time. For a classical spin, the 
main contribution comes from the derivative with respect 
to the component of the field that is perpendicular to the 
anisotropy axis || . The effect of the longitudinal field H z 
is expected to be small. The reason is that, classically, 
the energy scale over which r changes is determined by 
the very large anisotropy field Hk — 10 5 Oe. Note also 
that, classically, d 2 r/dH 2 \ < (see the inset of Fig- E} - 



For a molecular cluster, however, the magnetization 
reversal takes place by quantum tunneling. As in a two- 
level problem, the probability of tunneling is determined 
by the ratio between the tunnel splitting A and the bias 
£ generated by H z . For states with m — ±4 of Mni2 
we have A ~ 10 _1 K. As a result, r increases sharply 
between H z = and a few hundred Oe (<C Hk), as the 
inset of Fig. |2| shows Q. The large "positive" contribu- 
tion to X3 that we observe near lot = 1 can therefore 
be associated to the strong dependence of the tunneling 
probability on H z , which renders d 2 r/dH 2 \o > 0. 

In Fig. the experimental data are also compared 
with theoretical calculations that take into account the 
quantum nature of the molecular spins 0, ^| . The re- 
laxation time and susceptibility are calculated by solving 
a standard Pauli master equation, which includes tran- 
sition rates between energy eigenstates induced by the 
interaction with the phonon bath. Emulating the experi- 
mental conditions, we have computed \' an d x" f° r differ- 
ent applied fields \H\ > 50 Oe and fitted them to Eq. 0J. 
Despite the fact that the calculation does not take into 
account the precession of the spins, which produces the 
main classical contribution, it accounts remarkably well 
for the experimental results, including the change of sign 
with respect to the classical prediction. 

Before concluding, let us briefly consider the limit of 
large dc fields. Then the relaxation time approaches its 
classical limit, because the field blocks tunneling via all 
levels lying in energy below the classical anisotropy bar- 
rier. In order to check how this change of the relax- 
ation process modifies the nonlinear response, we have 
extracted, using the same method described above, \3 
from experimental data measured at \H\ > 1.5 x 10 3 Oe. 
The results, plotted in Fig. 01 resemble indeed the pre- 
dictions for classical spins. In particular, the inequality 
|%3 1 < Ix3t| is fulfilled in the whole frequency range. 

Interestingly, these high-field results agree well with 
calculations for overdamped classical spins (A > 1 in 
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FIG. 4: Real (•) and imaginary (o) nonlinear susceptibilities 
of Mni2 acetate obtained as described in the text [see Eq. Q] , 
from susceptibility data measured for \H\ > 1.5 x 10 3 Oe. 
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the upper panel of Fig. |3J). This is at first surprising 
since classically A w 2-it/ujlTq, where lot, — g^sHx/h 
is the Larmor frequency. For Mn 12 clusters, 2-k/oj l ~ 
4 x l(T 12 s while r ~ 1.3 x 10~ 8 s and 2.7 x 10- 8 s 
for Mni2Cl and Mni2 acetate, respectively. This gives 
A ~ 10~ 4 . For such a low damping the classical contri- 
bution is expected to dominate over the quantum one, 
whereas clearly the opposite is observed. 

We argue, however, that the situation might become 
different if the spins relax via quantum tunneling. Their 
precession around the anisotropy axis, as well as any 
other coherent dynamics, are then limited not by the de- 
cay towards lower lying energy states, but rather by the 
tunneling time r tun = h/A 17]. For times longer than 
Ttun the system evolves towards an incoherent mixture of 
delocalized (i.e., with (S z ) ~ 0) energy eigenstates. Near 
the top of the energy barrier (i.e., \m\ < 4), we have 
Ttun ^5 5 x 1CP 11 s. Replacing r tun for To in the above A 
gives A ~ 0.1, approximately equal for both compounds. 
Under these conditions, %3 is expected to depend weakly 
on the spin-lattice interaction, which agrees with the ob- 
served features. In addition, our results show that coher- 
ence is lost for times much shorter than to- 

In summary, we have investigated experimentally the 
nonlinear dynamical susceptibility of molecular nano- 
magnets and performed theoretical calculations at both 
the classical and quantum levels of description. We have 
found that Mni2 single-molecule magnets show a large 
quantum nonlinear magnetic response, which can be di- 
rectly linked to the tunneling mechanism by which their 
spins flip, and whose sign is opposite to the classical pre- 
dictions. This phenomenon should not be restricted to 
Mni2, but it could be observable as well for other meso- 
scopic magnetic systems (e.g., other magnetic clusters 
and nanoparticles) . The nonlinear magnetic suscepti- 
bility therefore offers a privileged probe to ascertain if 
tunneling plays a role in these mesoscopic materials, for 
which standard (in particular, linear) relaxation experi- 
ments do not always provide a definite answer [2(| . 
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